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Abstract 

In supersymmetric models with lepton-number violation, hence also R-parity vi- 
olation, it is easy to have realistic neutrino masses, but then leptogenesis becomes 
difficult to achieve. After explaining the general problems involved, we study the 
details of a model which escapes these constraints and generates a lepton asym- 
metry, which gets converted into the present observed baryon asymmetry of the 
Universe through the electroweak sphalerons. This model requires the presence of 
certain nonholomorphic R-parity violating terms. For completeness we also present 
the most general R-parity violating Lagrangian with soft nonholomorphic terms and 
study their consequences for the charged-scalar mass matrix. New contributions to 
neutrino masses in this scenario are discussed. 



PACS numbers: 12.60.Jv, 11.30.Fs, 98.80. Cq, 14.80.Ly 



1 Introduction 



The creation of a lepton asymmetry, i.e. leptogenesis |l|, ||, ||, which gets converted into the 
present observed baryon asymmetry of the Universe, is closely related to the mechanism by 
which neutrinos obtain mass. In general, all models of Majorana neutrino masses with the 
same low-energy particle content as that of the standard model are equivalent in the sense 
that they are all characterized by the same nonrenormalizable dimension-five operator [|J 
hr l ViVj(f) Q (j) . Different models of neutrino mass are merely different realizations || of this 
operator. They become distinguishable only at high energies, and since their interactions 
must violate lepton number, leptogenesis is a very natural possibility. For the canonical 
seesaw mechanism and the Higgs triplet model [0], leptogenesis does indeed occur nat- 
urally (l], 0. On the other hand, if neutrino masses are obtained radiatively ||, not only is 
leptogenesis difficult to achieve, the mechanism by which the former is accomplished leads 
naturally to the erasure of any primordial baryon asymmetry of the Universe M, EL ETUI. 



This is especially true in supersymmetric models of neutrino mass ITTJ with R-parity vio- 
lation. In a recent article fll2| , it was pointed out that leptogenesis is still possible in this 
case, provided that certain conditions regarding the R-parity violating terms are satisfied. 
Here we study this model in detail. 

In Section 2 we write down the superpotential of the lepton-number violating (but 
baryon-number conserving) extension of the supersymmetric standard model, together with 
all possible soft supersymmetric breaking terms, including the nonholomorphic terms [13|. 
In Section 3 we consider bilinear R-parity violation and how leptogenesis is related to neu- 
trino mass in this limited scenario. We find it to be negligible for realistic values of m v . 
In Section 4 we discuss how leptogenesis may occur without being constrained by neutrino 
mass in an expanded scenario. We assume negligible (enhanced) mixing between doublet 
(singlet) sleptons and charged Higgs bosons by allowing nonholomorphic soft supersym- 
metry breaking terms. In Section 5 we present the details of our calculations using the 
Boltzmann equations for obtaining the eventual lepton asymmetry. In Section 6, the com- 
plete charged- scalar mass matrix is displayed and analyzed. In Section 7, a new two-loop 
mechanism for neutrino mass is proposed. Finally in Section 8, there are some concluding 
remarks. 



2 Superpotential and Soft Supersymmetry Breaking 

In an unrestricted supersymmetric extension of the standard model of particle interactions, 
the chiral scalar superfields allow baryon-number violating terms which are not necessarily 
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suppressed. These dangerous terms are usually avoided by assuming a conserved discrete 
quantum number for each particle called R-parity, which is defined as 



R= (_ 1 )3B+L+2J j (1) 

where B is its baryon number, L its lepton number, and J its spin angular momentum. 
With this definition, the standard- model particles have R = +1 and their supersymmetric 
partners have R = — 1. We can list the three families of leptons and quarks of the standard 
model using the notation where all superfields are considered left-handed: 

Li = (ui, e iL ) ~ (1, 2, —1/2), e\~ (1,1,1), (2) 
Qi = (uudi) ~ (3,2,1/6), (3) 
<~ (3*, 1,-2/3), d?~(3*, 1,1/3), (4) 

where i is the family index, and the two Higgs doublets are given by 

#i = (/i?,/ir)~(l,2,-l/2), (5) 
H 2 = (h+,h° 2 ) -(1,2,1/2), (6) 

where the SU (3)c x 5/7 (2) ^ x /7(l)y content of each superfield is also indicated. If R-parity 
is conserved, the superpotential is restricted to have only the terms 

w = iie x r 2 + r.u.i-y] + .n'lhQ,'!) + f^Q^. (7) 

In this case, both baryon and lepton numbers are conserved. However, to forbid proton 
decay, it is sufficient to conserve either baryon number or lepton number (because the final 
state of the proton decay must contain a lepton or antilepton). If only baryon number 
or only lepton number is violated (thus R-parity is also violated), the conservation of the 
other quantum number is enough to satisfy all present experimental constraints. This has 
motivated numerous studies of R-parity violating models. 

If R-parity is violated, the superpotential contains the additional terms 

W = HiLiH 2 + XykLiLjel + X'^LiQjdl + Aj fe <d^. (8) 

We cannot have all of these terms because then the proton will decay very quickly. We 
may choose only the lepton-number violating terms or only the baryon-number violating 
terms. Following the overwhelming choice of many others, we consider only the former 
case and set \"- k = 0. The remaining terms may now induce nonzero neutrino masses, 



either from mixing with the neutralino mass matrix, or in one- loop order O]. Although 
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these terms are allowed, we do not know if they originate from any fundamental theory, so 
the couplings are considered free parameters, constrained only by experiment. 

Other free parameters exist in the minimal supersymmetric standard model (MSSM), 
i.e. the soft supersymmetry breaking terms, which do not introduce quadratic divergences 
to the renormalized theory. Usually only the holomorphic terms are considered which come 
from the chiral superpotential interacting with gravity, together with the gaugino masses. 
The most general such Lagrangian conserving R-parity is: 

C soft = -^(M^L" - e^M^ 

-~dt{M 2 d ) l3 ~d) - M 2 Hl H?H« - M 2 H2 H?H« - e ab (BH*H b + h.c.) 
Sab^Aef^HfL^ + {A u f u ) ij H b 2 Q'm + (A d f d ) ij HiQ b i d c j + h.c. ^ 

-- (M 3 gg + M 2 WW + M X BB + h.c.^j . (9) 

If R-parity is violated, more soft terms may be present, i.e. 

Cf oft = -e ab {B[LtH b 2 + A% k LfL b e c k + Af jk QfL b d c k ) - A% k u\d]d\ + h.c. . (10) 

We follow the convention that the coupling constants of all the R-parity conserving soft 
terms are denoted without a prime, while the R-parity violating terms are denoted with a 
prime. 

Since the soft terms may originate from gravity couplings, there is no clear reason 
for them to come only from the renormalizable chiral superpotential. They may also 
originate from nonrenormalizable terms, which can be functions of both left and right chiral 



superfields. Such nonholomorphic terms allow new supersymmetry-breaking soft terms |T3] 



in the Lagrangian. The most general set of nonholomorphic soft terms conserving R-parity 
is: 

C?o% = -N^HTlm - NtFTQfd) - ^BTQIu) + h.c. . (11) 
Similarly, nonholomorphic soft terms breaking R-parity are: 

CZS? = -N'?H^*L* - N- e H%*H"et - N-j k Lf*Q"u c k 

IS ~c~c~ic* TKT'd ^ r\o-r\bjc* 



~N i j k u c i c C jd c k - N% k e ab Q?Q°d c k * + h.c. (12) 

In our convention, all "N" constants are for nonholomorphic terms. Since we assume lepton- 
number violation but baryon-number conservation, it implies \'^ k = A'^ k = N'^ k = 0. 
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3 Bilinear R-Parity Violation 



Before coming to the explicit model in the next section, we first look at the possibility 
of generating a lepton asymmetry from bilinear R-parity violating terms. This case is a 
good illustration of the general problems involved. One immediate consequence of the 
violation of lepton number through the bilinear terms is the mixing of the neutrinos with 
the neutralinos. In the MSSM there are four neutralinos, the U(l) gaugino (B), the 
SU(2) gaugino (W 3 ), and the two Higgsinos (h®, h®)- When lepton number is violated 
through the R-parity violating terms, it is possible to assign zero lepton number to what 
we usually regard as lepton superfields [14], |15| . Suppose only the r neutrino mixes with 
the neutralinos, then both r and v T may be assigned an effective vanishing lepton number 
while the other leptons remain leptons. However, in the general three-family case, all three 
neutrinos may mix with the neutralinos, and the scalar partners of the neutrinos [pi) may 
all acquire nonzero vacuum expectation values (VEVs). 

In the most general case the neutralino mass matrix with all the seven fields in the 



basis 


B 


w 3 , hi 




V%y v 3 


is given 


by 










Mi 





-sr z vi 


sr z v 2 


-sr z v Vl 


-sr z v U2 


-sr z v V3 









M 2 


cr z v\ 


-cr z v 2 


cr z v Ul 


cr z v U2 


cr z v U3 






-sr z vi 


cr z vi 





~M 











M 




sr z v 2 


-cr z v 2 


-V 







— A*2 


-Ms 






-sr z v Ul 


cr z v Ul 





-Mi 















-sr z v U2 


cr z v U2 




















-sr z v vz 


cr z v V3 





— A*3 












(13) 



where s = sxn9w, c = cos9w, r z = M z /v, and v\, v 2 , v Vi are the VEVs of hi, h 2 , and 



(246 GeV) 2 and v 2 v 



v 2 + v 2 + v 2 . 



We also define 



respectively, with v\ + v 2 + v 2 = v 
tan/5 = v 2 /{y 2 + vl) 1 / 2 . 

To understand the structure of the above 7x7 mass matrix, let us assume that /x is 
the dominant term, then h\ 2 form a heavy Dirac particle of mass \x which mixes very little 
with the other physical fields. Removing these heavy fields will then give us the reduced 
5x5 matrix in the basis [B, W3, V\, v 2 , v 3 ): 



M 



Mi - s 2 5 


scS 


-sei 


-se 2 


-se 3 


scS 


M 2 - c 2 5 


cei 


ce 2 


ce 3 


-se 1 


cei 











se 2 


ce 2 











-se 3 


ce 3 












(14) 
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where 



r2 v lV2 l Mf sin 2(3 I v 2 



5 = 2M^-= Z ' \ (15) 
v 2 fi n y v 1 cos 2 p 

Mz ( »i \ c\a\ 
e t = —[v^-- Vi y (16) 

From the above, only the combination v\ = [e\V\ + e 2 ^2 + £3^3) /e, with e 2 = e 2 + e| + e§, 
mixes with the gauginos. This state will have an effective vanishing lepton number and 
the other two orthogonal combinations decouple from the neutralino mass matrix. In this 
case, only the eigenstate 

* = « + (17) 



gets a seesaw mass, i.e. 



s 2 c 2 



^^{k + k)' (18) 

whereas the other two neutrinos remain massless. They may get masses through one-loop 
radiative corrections from the usual trilinear R-parity violating terms which we have not 
yet considered. 

The two gauginos mix with the neutrino v% and form mass eigenstates given by 

~ ~ scS ~ sc 

B ' = B + i^w^-w^ < 19 > 

~ , ~ sc5 ce , 

Wl = W 3 -— —B + —ui. 20 

The physical states B' and W' z now contain v\. This gives the main feature of R-parity 
violation, which is the decay of the lightest neutralino. By virtue of their v x components, 
both neutralinos will now decay into a lepton or an antilepton and a weak gauge boson, 
such as W' z — > l~W + and l + W~ , thus violating lepton number. Since the mixing of the 
neutralinos may also have CP violation through the complex gaugino masses (thus making 
5 complex), a lepton asymmetry may be generated from these decays. However, the amount 
of asymmetry thus generated is several orders of magnitude too small because it has to be 
much less than (e/Mi^) 2 , which is of order m u > /Mi^, i.e. < 5 x 10~ 13 if r?v < 0.05 eV and 
M 12 > 100 GeV. In addition, the out-of-equilibrium condition on the decay width of the 
lightest neutralino imposes an upper bound on (e/Mi^) 2 which is independent of m^', and 
that also results in an asymmetry very much less than 10~ 10 . 

We now consider the R-parity violating trilinear couplings, i.e. A and A' of Eq. (8). Since 
the particles involved should have masses at most equal to the supersymmetry breaking 
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scale, i.e. a few TeV, their L violation together with the B + L violation by sphalerons 
[16| would erase any primordial B or L asymmetry of the Universe 0. To avoid such a 
possibility, we may reduce A and A' to less than about 1CU 7 , but a typical minimum value 
of 10~ 4 is required for realistic neutrino masses in one-loop order ||11|| . Hence it appears 
that the MSSM with R-parity violation is not only unsuitable for leptogenesis, it is also a 
destroyer of any lepton or baryon asymmetry which may have been created by some other 
means before the electroweak phase transition. 



4 Leptogenesis from Neutralino Decay 

From the discussion of the previous section we observe that for a leptogenesis mechanism 
to be successful in the MSSM with R-parity violation, two requirements have to be fulfilled. 
First we must use lepton-number violating terms which are not constrained by neutrino 
masses. Second we must satisfy the out-of-equilibrium condition for the lightest neutralino 
in such a way that the asymmetry is not automatically suppressed. More explicitly, we will 
consider the possibility that the heavier neutralino does not satisfy the out-of-equilibrium 
condition and decays very quickly, but the lighter neutralino decays very slowly and satisfies 
the out-of-equilibrium condition. Since the asymmetry comes from the interference of the 
one-loop CP violating contribution of the heavier neutralino, it is then unsuppressed. We 
will demonstrate explicitly in the following how this scenario may be realized. 

We assume first that Mi > M 2 , so that the bino B is heavier than the wino W3. While 
the former couples to both e iL e iL and e?e?, the latter couples only to e iL e iL , because the 
el are singlets under SU{2) L . Since R-parity is violated, one combination of the e iL and 
another of the e~l mix with the charged Higgs boson of the supersymmetric standard model: 

= hf cos (3 + hf sin (3. Let us denote them by II and l c respectively. Their corresponding 
leptons are of course II and l c . Hence both B' and W3 may decay into Z =F /i ± . 

We assume next that the II mixing with h~ is negligible, so that the only relevant 
coupling is that of B to l c h + . Hence W% decay (into l^h^) is suppressed because it may 
only do so through the small component of B that it contains, assuming of course that all 
charged sleptons are heavier than B or W3. 

With this choice that the heavier neutralino B' decays quickly and the lighter neutralino 
W% decays much more slowly, we now envisage the following leptogenesis scenario. At 
temperatures well above T = Msusy, there are fast lepton-number and R-parity violating 
interactions, which will wash out any L or B asymmetry of the Universe in the presence of 
sphalerons. This will be the case even at temperatures around Mi, when W$ interactions 
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violate Li as well as B—3Li for i = e, /i, r 0. We assume here that all other supersymmetric 
particles are heavier than the neutralinos, so that at temperatures below Mi we need only 
consider the interactions of B' and W%. In Figure 1 we show the lepton-number violating 
processes (a) B' <-> l%h T , where we have adopted the more conventional notation of an 
outgoing Ir in place of an incoming l c . These processes are certainly still fast and there can 
be no L asymmetry. At temperatures far below the mass of the heavier neutralino, the & 
interactions are suppressed and we need only consider those of W%. With our assumptions, 
the lepton-number violating processes (b) W' z <-> l%h T are slow and will satisfy the out- 
of-equilibrium condition for generating a lepton asymmetry of the Universe. Specifically, 
it comes from the interference of this tree-level diagram with the one-loop (c) self-energy 
and (d) vertex diagrams. Since the unsuppressed lepton-number violating couplings of B' 
are involved, a realistic lepton asymmetry may be generated. It is then converted by the 
still active sphalerons into the present observed baryon asymmetry of the Universe. In 
this scenario the mass of also has to be small enough so that the scattering processes 
mediated by the heavier B' are negligible at temperature below M 2 when the asymmetry 
is produced. 

We start with the well-known interaction of B with I and Ir given by flT 

eV2 



cos 9 



IV 



I ( I BI H - //, . 



(21) 



We then allow Ir to mix with h , and B to mix with W3, so that the interaction of the 
physical state W$ of Eq. (|20| ) with I and is given by 



sc£5r \ I 



Mi - M 2 J \cos9 



(22) 



where £ represents the Ir — h~ mixing and is assumed real, but the parameter 5 of Eq. (15) 
is complex. We have also inserted a correction factor r = (1 + M 2 /fi sin 2(3) /{I — Mf/yU 2 ) 
for finite values of M 2 / fi. The origin of a nontrivial CP phase in the above is from the 2x2 
Majorana mass matrix spanning B and W%, with complex Mi and M 2 . It is independent 
of the phase of \x and contributes negligibly to the neutron electric dipole moment because 
the magnitude of 5 is very small. [Note that the usual assumption of CP violation in 
supersymmetric models is that Mi and M 2 have a common phase, in which case the phase 
of 5 would be equal to the phase of fj,.] 
The decay width of the bino is then 

r* = r(& - + r(& -> ™+) = (23) 
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while that of the wino is 



47T \Mx-M2 



2\2 



(24) 



Using Eqs. (pip and (|22|), we calculate the interference between the tree- level and self- 
energy + vertex diagrams of Figure 1 and obtain the following asymmetry from the decay 



of W^. 



a£ 2 im<5 2 ( 1 _ _m|_ y x^(x) g 



2cos 2 # w 



M 2 T 



W3. 



(1-x) 



where a; = M 2 - ,/M|, and 

' B' 



<7(x) = 1 + 



2(1 - x) 



X 



1 + x 



x 



ln(l + x) -1 



(26) 



If the W% interactions satisfy the out-of-equilibrium condition, then a lepton asymmetry 
may be generated from the above decay asymmetry. Note that in the above expression for 
e, the parameter 5 appears only in the combination Im<5 2 /|5| 2 , which may be of order one. 
If the absolute value of 5 is small, it slows down the decay rate of and a departure from 
equilibrium may be achieved without affecting the amount of decay asymmetry generated 
in the process. 



At the time this lepton asymmetry is generated, if the sphaleron interactions [16| are 
still in equilibrium, they will convert it into a baryon asymmetry of the Universe 0. If the 
electroweak phase transition is strongly first-order, the sphaleron interactions freeze out 
at the critical temperature. Lattice simulations suggest that for a Higgs mass of around 
ran ~ 70 GeV, the critical temperature is around T c m 150 GeV [18]. Higher values of ran 
will increase the critical temperature, but the increase is slower than linear. For example, 
for van ~ 150 GeV, the critical temperature could go up to T c w 250 GeV. 

For a second-order or weakly first-order phase transitionF], the sphaleron interactions 
freeze out at a temperature lower than the critical temperature. After the electroweak 
phase transition (T < T c ), the sphaleron transition rate is given by [p0| 



T sph (T) = (2.2 x 10 4 k) 



[2 M W (T)]< 
[4tt a w T] 3 



IT 



(27) 



1 Note that we do not require the electroweak phase transition to be first-order for satisfying the 
out-of-equilibrium condition. See for example Ref. fj~9|j , 
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where ^ 

M w (T) = \g 2 (v(T)}= 1 -g 2 (v(T = 0)} (l - ^ , (28) 

and the free energy of the sphaleron is E sph (T) « (2 M W /(T)/a(W / )) B(m,h/M w ), with 
-B(O) = 1.52, B(oo) = 2.72 and k = e -3 ' 6 ||21|| . In this case, the sphalerons freeze out at a 
temperature T out which is the temperature at which their interaction strength equals the 
expansion rate of the universe, 

T sph (T out ) = H(T out ) = 1.7V^|k (29) 

Mpi 

For a critical temperature of about T c ~ 250 GeV, the freeze-out temperature comes out 
to be around T out ~ 200 GeV. 

These discussions indicate that as long as the lepton asymmetry is generated at a 
temperature above, say 200 GeV, it will be converted to a baryon asymmetry of the Uni- 
verse. Since the sphaleron interactions grow exponentially fast, they can convert a lepton 
asymmetry to a baryon asymmetry [a (B — L) asymmetry to be precise] by the time 
the temperature drops by only a few GeV. In the next section we discuss how the decay 
asymmetry of the neutralinos becomes a lepton asymmetry of the Universe. 



5 Boltzmann Equations 



We now solve the Boltzmann equations [[24J] to estimate the amount of lepton asymmetry 
created after the decays of the neutralinos. When the decay of the satisfies the out- 
of-equilibrium condition, i.e. when the decay rate is slower than the expansion rate of 
the Universe, the generated asymmetry is of the order of the decay asymmetry given in 
Eq. (P5j). This argument could replace the details of solving the Boltzmann equations 
for an order-of-magnitude estimate of the asymmetry in many scenarios. However, in the 
present case there are other constraints and depleting factors, and we need to solve the 
Boltzmann equations explicitly for a reliable estimate. 

If the decay rate is much less than the expansion rate of the Universe, the gener- 
ated lepton asymmetry is the same as the decay asymmetry. In other words, the out-of- 
equilibrium condition reads 

K * -wkr)^ (30) 

where H(T) is the Hubble constant at the temperature T and is given by 



H[T) = /2 . (31) 

" 45 M Planck ' 1 ; 
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with the number of massless degrees of freedom which we take equal to 106.75 and 
Mpi ~ 10 18 GeV is the Planck scale. If this condition is satisfied, the lepton asymmetry 
is given by til = ni — nj ~ e/g*. But in practice, when K <C 1, there is no time for the 
asymmetry to grow to its maximum value before the sphaleron transitions are over. So we 
need to study the case K ~ 1. Furthermore, a reasonable amount of asymmetry cannot be 
obtained unless the inverse decay and the scattering from bino exchange have rates lower 
than the expansion rate of the Universe. All these effects result in the further diminution 
of the lepton asymmetry and we need to solve the Boltzmann equations to take care of 
them properly. 

At temperatures T < M 2 , the decays of W% given in Eq. (|24| ) start generating an 
asymmetry. At this time there are important damping contributions coming from the 
inverse decays of W% and B' as well as the scattering processes + h T — ► B' — ► l T + . 
As we will see, the last two processes are especially important because B' tends to remain in 
equilibrium and its presence washes out the created lepton asymmetry from the decays. 
The reason is that its interactions are strong enough so that the Boltzmann exponential 
suppression of its number density may not be sufficient to compensate its large inverse 
decay and scattering cross sections. The effect of the scattering + h T — > W' z — > l T + 
is on the other hand negligible because it is suppressed by a factor of [(sc5r) / \M 1 — M 2 )} 2 
with respect to the scattering I + h T — > B' — ► l T + h . Neglecting this term and defining 
the variable z = M^/T, the Boltzmann equations are then: 



^3 = eq ~ ""3 i| (on) 

dz ^ MsH{M w ,)\X e l l 



dX 



L 



eq 

V 



dz 'wi s H(Mi 



3 

X WL \ 1 X 



Xfl J 2 X. 



L 



sH(M B ,)\M w ,. 



(33) 



where we have defined the number densities per comoving volume X{ = rii/s in terms of 
the number densities of particles "i" and 

2vr 2 

S = 9*^T 3 (34) 

is the entropy density. The equilibrium distributions of the number densities are given by 
the Maxwell-Boltzmann statistics: 



M~ 
w. 



w' 



9W -^±TK 2 (M W ,/T) } (35) 



3 2vr 2 
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M 2 ~ 

n B> = 9b>^-TK 2 (M § ,/T), (36) 
n 7 = (37) 

where = 1, ^ = 1, and g 7 = 2 are the numbers of degrees of freedom of PVg, B', and 
the photon respectively. 

The quantities 7^, and 7J? are the reaction densities for the decays and inverse decays 



of W£ and B'\ 



K\(M^r,/T) 

<y eq — n eq 3 V - , 

^w> ~ n wiK 2 (M w ,/T) W ^ W 



eq _ n e,W) r . m 

~ n B'K 2 (M § ,/T) LB " [6J) 

Ki and K 2 being the usual modified Bessel functions. The reaction density for the scat- 
tering is given by 

Islatt. = 7T— 4 / dsa §l ( S ) yfsK^yfs/T), (40) 

^ J(m h +mi) 2 

where <7g, is the reduced cross section and is given by 2[s — (m^ + m/) 2 ][s — (rrih — 
m i) 2 } (J B'/ s ~ 2sa B ,. The cross section a B , does not contain the contribution of the 
on-mass-shell bino (which is already taken into account in the decay and inverse decay 
terms). This is achieved by replacing the usual propagator l/(s — m? + iTm) with the 
off-mass-shell propagator |22|. ^3| : 



s — m 2 



D' 1 = — — . (41) 

(s — m 2 ) 2 + T 2 m 2 

The cross section a B , contains the s- and t-channel contributions together with their inter- 
ference terms and is given by 



a B , = a^h* -> & -> l T h+ 
1 / e 2 e \2 



87TS 2 \COS 2 8 w J B 



s 2 As 2s / s + Mlx / s 



(42) 



In Eq. PD]) the integral is dominated by the s-channel contribution in the resonance region 
and to a very good approximation, %catt. reduces to 



^ = 512^ Tj™'/ T >^ (43) 
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From Eqs. fl3"3|) to (]43|), we find the Boltzmann equations, i.e. (|32| ) and (0), to be given 
by 

dz ~ w *K 2 {z)\ w * W3 



dX L r ^ K x (z 



1 X - 



A 



2K 2 (zM B ,/M w ,) X y " X 7 sT^, 



(44) 



with 



K & = r §l /H(Ma,), (45) 

which gives the strength of lepton-number violation in the decays of B' . 

If we now ignore the inverse decay and scattering processes, we can simplify the problem 
by requiring the out-of-equilibrium condition to be 

K w , < 1. (46) 

With the terms proportional to Kg, this condition is necessary but not sufficient. In the 
present scenario for a large asymmetry we also require Kg, > 1. Indeed, K B , is larger 
than K w , by a factor R K = K^/K^ ~ \(sc5r)/(M l - M 2 )]~ 2 (M W ,/M B ,), which is larger 
than one by several orders of magnitude. Therefore B' remains in equilibrium and the 
B' damping terms, due to its inverse decay and scattering, dominate over the W% inverse 
decay damping term as long as the Boltzmann suppression factor in the B' equilibrium 
distribution has not compensated the large value of Rk- For example with the set of 
parameters M w , = 2 TeV, M& = 3 TeV, sin 2(3 = 0.5, ( = 2x lCT 3 , \l = 5 TeV used in 
Ref. [12], we obtain K^, = 0.63 and Kg, = 7.8 x 10 5 , and the B' damping terms dominate 



over those of W% as long as the temperature is above ~ 65 GeV (the former differs from 
the latter by a factor ~ R K {M B ,/M w ,fe~ {M& ~ M ^ )/T ). In this case the inverse decay and 
scattering of & cause a considerable wash-out of the asymmetry because W% has mostly 
decayed away already at temperatures well above ~ 65 GeV. This is illustrated in Figure 
2 showing the effects of various terms in the Boltzmann equations. 

To avoid this wash-out, the value of Afg, has to be larger in order that the B' number 
density is further suppressed at temperatures below M-^, when the asymmetry is produced. 
Varying the parameters of these B' damping terms, it appears difficult to induce a suffi- 
ciently large asymmetry of order 10~ 10 for M B , below 4 TeV. Two typical situations for 
which a large asymmetry is produced are for example: 

M B , = 6 TeV, M w , = 3.5 TeV, f = 5 x 10~ 3 , 
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H = 10 TeV, sin 2(3 



0.10, m h = 200 GeV, 



(47) 



and 



M^, = 5TeV, M w , 
/i = 7.5 TeV, sin 2(3 



2 TeV, f = 5 x 10~ 3 , 
0.05, m h = 200 GeV, 



(48) 



for which we have K w , = 0.02, K' & = 2.4 x 10 6 and K w , = 0.02, K'~, = 2.9 x 10 6 
respectively. At T = Mz the leptonic asymmetry produced is = 1.0 x 10~ 10 with the 
parameters of Eqs. ([I7|) and Xi = 1.2 x 10~ 10 with those of Eq. (f4*8"D. Figures 3 and 4 show 
the evolution of the asymmetry in these two cases. As can be seen from these figures, the 
damping effects of the inverse decay of and of the scattering are smallf]. The damping 
effect from the inverse decay of B' is however not small and reduces the asymmetry by a 
factor of 2 to 4 by washing out all the asymmetry produced above T ~ 300 — 400 GeV. 

A large asymmetry of order 10~ 10 is produced provided M^, is of the order 4 TeV 
or more. A low value of sin 2(3 below ~ 0.30 is generally necessary. Values of £ around 
3 — 5 x 10~ 3 , of /i around 5 — 10 TeV, and of M-^, from 1 TeV to 2M^,/3 are also preferred. 

6 Charged Scalar Mass Matrix 

The mechanism we propose for leptogenesis requires the decay of B' to be fast, while that 
of W% is very slow. This is achieved by requiring l^h T to be the main decay mode and 
l L h T to be negligible. Hence Ir must mix with h~ readily, so that B' could decay directly, 
but could decay only through its small B component. We now consider the charged 
scalar mass matrix which determines this mixing. As shown in the following, our present 
scenario requires one more new ingredient, i.e. the presence of nonholomorphic soft terms. 

The value of the tr — h + mixing parameter £ is governed by the charged scalar mass 
matrix which follows from the quadratic terms in the Lagrangian: 



with $ = [h^*,h+,e-*,e'r} T . In the case where all N constants are put to zero and for 
parameters satisfying the various constraints from the Boltzmann equations (see previous 
section), it is difficult to generate a sufficiently large value of the l c — h~ mixing param- 
eter £ (see Appendix A). Thus we will ignore the /ij terms and the associated vaccuum 

2 Note that in the Boltzmann equations we neglected the clamping contributions of the scatterings 
V^P 1 — > /i ± /i ± mediated by a neutralino in the t-channel. Their effect is negligible within the ranges of the 
parameters we consider. 



£ 3 -&M 2 S ±$, 



(49) 
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expectation values v Ui in the following. To induce a large e c — h mixing, we introduce the 
nonholomorphic terms of Eqs. (|TT|) and (12). The mass matrix is then given by 



M 2 



s= 



M 2 



-hi 



M 2 



i- h T 



Ml 
Ml 

tin 

M 2 
M 2 



-h+ 



L J L 2 



Ml 
Ml 

tin 

M 2 
Ml 



e,'. — e* 



Ml 
Ml 



Ml 
M 



JL 
2 



(50) 



where 



M 2 ^ 


-hi 




-h+ 


M 2 

h i - 








M\- 


-n 


M 2 

h 2 - 


-h+ 






Ml- 


~% 


M 2 . 




M 2 . 


-n 


Ml* 

3 





9 2 2 

1 V2 

9 2 



B 



V2 
Vl 

B. 



n: 



IB 



g 2 2 u v i 

— V7 — B — , 

4 1 v 2 ' 

-B' 



V2 



11 1 

mh - 9 -^{v\- v 2 2 )5 l3 + ~r ki r kj vi 



(51) 

(52) 
(53) 

(54) 

(55) 
(56) 

(57) 

(58) 
(59) 

(60) 



In the following we assume that N- B and B[ are negligible so that they do not induce 
a large mixing of the left-handed charged sleptons with the charged Higgs boson. We are 
then left with the N- e and terms. Going to the basis of the physical charged Higgs 
boson h + and of the Goldstone boson G + , the latter decouples and in the basis [/i H 
we obtain the mass matrix 



> e i L , e i\ 



( 



M 2 



m 2 w -2 



sin 2/3 







Ml 



M 



3L~ 
2 



e.—e. 
3 *L 



T2 N " v ^ 

M%_~ e c 



(61) 
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We observe that only one combination of the charged right-handed sleptons mixes with 
the charged Higgs boson: 

N' e P c -I- NL e Z? -I- NL e £c 

(62) 



N ,e 

with (N' e ) 2 = (N[ e ) 2 + {N2) 2 + (N3) 2 . In the Lagrangian the mass term which couples 
the charged Higgs boson with the sleptons reduces therefore to the following single term: 



C 3 — l -vN' e h + l c + h.c. 
y/2 



(63) 



In the case of one family, we obtain to a very good approximation: 

£ = 



T2 N ' ev 



m 



h+ 



M? 

e L -e* 



m 2 ,-M? _ 
h+ e L~ e L 



Mi-* 



III 



h+ 



(64) 



(65) 



where £' is the h + — mixing and 



B 



ni 



h+ 



m 



sin 2(3 



(66) 



As required in section 3, the mixing £' has to be much smaller than the mixing £ in order 
to avoid having the transition W' z — > e^ej — > ^ ± ej. This requires 



2s 2 5r 



Mi - M 2 



(67) 



which for the parameters of Eq. Q4"7D implies that f < 7x 10~ 5 £. Our mechanism requires 
also that 771^+ < M^ and that the mass of any charged slepton is larger than M^. We 
then have 



mr 



M? 

e L ~e L 



m e \i tan (3 + A e m e H — —N e v 2 

V2 



< Mi 



2s 2 5r 



Mi - M 2 



(68) 



(69) 



which for the set of parameters given in Eq. (|47|) and for M ( 



M|, gives iV /e ~ 1 TeV 
and |m e //tan/3 + A e m e + ^A^" e t> 2 | < (50GeV) 2 . The latter condition requires m e < 0.02 
GeV, \A e m e \ < (50 GeV) 2 , and |iV e | < 15 GeV, or a cancellation of the three terms together. 
In the case where there is no cancellation between these terms the charged slepton which 
mixs with the charged Higgs boson must have predominantly an electron or /x flavor. In 
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the case of a cancellation between these terms, all flavors are possible. For other sets of 
parameters which lead to large asymmetries, these numerical bounds can be relaxed easily 
by a factor of 2 to 4. In the more general case of three families, similar constraints and 
relations are obtained. 



7 Two-Loop Neutrino Mass 

It is interesting to note that in addition to inducing a lepton asymmetry, the nonholomor- 
phic terms N- e could also generate a neutrino mass. Since lepton number is violated at most 
by one unit in each term, the neutrino mass should include at least two lepton-violating 
vertices in a loop diagram. 

There exists a one-loop diagram (Fig. 5), which contributes to the sneutrino "Majorana" 
mass. In general, the sneutrinos could have diagonal lepton number conserving masses, i.e. 



v*v. But there can be also lepton-number violating mass terms, i.e. vv p5| , Pq , |27| . In the 
present model, the nonholomorphic terms give rise to a lepton-number violating sneutrino- 
antisneutrino mixing term, i.e. C 3 —^Srn^.MiPj + h.c. In the case of one family, we 

get 

5m l ~ ^— m i- ( 7 °) 
This lepton-number violating sneutrino mass can then induce a Majorana neutrino mass 
2|. 

In the present case the one-loop diagram of Fig. 6 gives a neutrino mass 



1 e 2 2 
m v ~ — — r . o „ orrirMw, 



M? — M - , — M 2 - , In ( M? IM 2 ~ , 



32^^^ " -3 (Ml-M 2 ,,) 2 



256vr 4 sin 2 ^ v 2 w * ( M 2 _ M 2 

In the case where the lepton I which mixes with h + is essentially r, we get 

m ~ ~ 256^^6^^^ (m? -Ml,) 2 ' 

which has the correct order of magnitude. For example with the parameters of Eq. ( |47D 
and taking My ~ Mg, we get m„ T ~ 0.1 eV. The value of ^ we need for having the right 
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order of magnitude for the asymmetry is therefore also the one we need to have a neutrino 
mass, in agreement with the present data on atmospheric neutrinos. In Eq. (|7l|) , the 
factor 8ml appears because there is GIM (Glashow-Iliopoulos-Maiani) suppression from 
summing over all possible neutral slepton eigenstates in the loop. In Fig. 5 the two- 
point functions of the form f(m 2 , m' 2 , p) = {i/ir 2 ) J d 4 k(k 2 — m 2 )~ 1 ((/c + p)' 



m 



/2\ 



have 



been (roughly) approximated by ~ 1 while in Fig. 6 the two-point functions have been 
calculated explicitly (as required by the fact that for these diagrams, a GIM suppression 
mechanism is operative). This can also be understood from another point of view. Since the 
diagonal terms of the sneutrino mass come from the lepton-number conserving interactions, 
they should not contribute to the Majorana mass of a neutrino. Only the lepton-number 
violating sneutrino mass, which is the mass-squared difference, should contribute to the 
Majorana neutrino mass. This makes the neutrino mass proportional to the mass-squared 
difference after GIM cancellation. Note that combining both one-loop diagrams, a two-loop 
diagram is obtained which is similar to the diagram proposed in Ref. |28j with a different 
lepton-number violating soft term. 

In the case of three families, the induced neutrino mass terms involving v e are suppressed 
by the small value of the electron mass with respect to the /i or r mass. Therefore, unless 
N' e e is much larger than N'^ and Njf, u e essentially decouples and acquires a very small 
mass; we get m Ue ~ 10 -8 eV or less. In this case, the mass matrix of the sneutrinos in the 



IM — T sector is of the form: 



C 3 -IqImI*,, 



(73) 



with 



r, v\, v\) T and 



Mr 



( M 2 



5m 2 - 



Ml 
8ml 



\ 8ml 8ml 

> "lit ^TT 



8ml 
8ml 

Mr 





8ml \ 
8ml 


Ml ) 



(74) 



where for simplicity we have assumed in Eq. @ a diagonal matrix M\ = diag{M\ , M 2 t ] 
and with 

,2f2 N' e N' e 



8m-.. ~ 



btt 2 v 2 mumh {N'*y 



(75) 



From the mass matrix of Eq. ([74]), the diagrams of Fig. 6 induce then the following neutrino 
mass term 

C 3 -hj/lM^ v , (76) 
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with ty v = fa, ujt, u T , v\) T and 



where 



M v ~ A 



( A M _ M A„_ T \ 

A^ A^„ T 

A M _ T A r _ r 

\ A M _ T A T _ T / 



(77) 



A 



and 



M w ,(f(Mt + Sm^Ml 0) - /(M£ - M^,, 0) 



2<W,M^, 



M 2 ~ 



Mi 



Mr 



M W , f{Ml , M? , 0) - f(Mi , M^ , 0) 



^rn%M w ,{Ml - M^)" 1 ^ - M^,)" 1 ^ - M^)" 1 

M? 



Mf M| , In 



Mr 



M 



+ Ml T M? In 



M 2 M 2 In 



A 



M 



M] 



64?r 2 sin 2 6 W ' 

This matrix can lead easily to a maximal mixing between the \x and r neutrinos, 
will be the case in particular if A^.^ ~ A T _ r which implies 



m T N' T e 



Mf ' 



,(79) 

(80) 
This 

(81) 



M| , the determinant of 



In addition it can be seen easily that in the limit where M 
the neutrino mass matrix vanishes, leading to a large hierarchy of masses (as required by 
atmospheric and solar neutrino experiments, taking into account the fact that the mass of 
v e is below 10~ 8 eV in the present scenario). For example with the parameters of Eq. ( f47|) 



and taking in addition N'* = UN' T e , M Lm ~ = 6 TeV, and M Lr ~ 7.5 TeV, we obtain 
one neutrino with a mass ~ 10~ 3 eV and one with a mass ~ 10~ 5 eV in addition to the 
electron neutrino with a mass below 10~ 8 eV. In this case the mixing between the \x and r 
flavors is large (sin 2a = 0.99) while that of the electron flavor with the two other flavors 
is very much suppressed. Note that the values of the lepton-number violating mass terms 
5m|. . induced by Fig. 5 are several orders of magnitude below the phenomenological bounds 
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5m e < 350 MeV, 5m^ < 50 GeV, and 5m T < 450 GeV obtained for M SUSY ~ 1 TeV in 



Ref. [27 



In summary, from the above qualitative estimate, we observe that realistic neutrino 
masses could be accomodated easily in the present scenario, in agreement with atmospheric 
and solar neutrino experiments. A large mixing and a hierarchy of neutrino masses appear 
rather naturally. A more quantitative estimate would require an explicit calculation of the 
two-loop integrals involved, but since there are still many free parameters, it will not add 
much to our understanding in any case. 



8 Conclusion 

We have studied a model of leptogenesis in a R-parity violating supersymmetric model. 
The lightest neutralino W% is assumed to be mostly the 577(2) gaugino but its decay into 
/ ± /i T is suppressed because the required l^— h + mixing is negligible. On the other hand, W% 
has a small component of B, the U(l) gaugino, which decays readily because the required 
I Ft ; — h + mixing is of order 10 -3 from the presence of nonholomorphic R-parity violating 
soft terms in the Lagrangian. The decay asymmetry of W% is then evolved into a lepton 
asymmetry of the Universe by solving the Boltzmann equations in detail numerically. We 
demonstrate how each term in the equations affects the eventual outcome of the proposed 
scenario. The charged scalar mass matrix and the neutralino sector are discussed in detail. 
A realistic scenario of radiative neutrino mass generation in two loops is presented, which 
originates from the same lepton-number violating nonholomorphic terms. 
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A Complete Charged Scalar Mass Matrix 



From Eqs. (0) to (|12D, neglecting small terms of order (A^) 2 and Xijk^'imn, the com- 
plete charged scalar mass matrix is given by Eq. ([5(]) with (see also Refs. [15|, |2£| for the 
holomorphic part): 



M 2 h - _ c 

h, —e c 
l % 



Ml 



M 



h n — e* 



Ml 



M\ 



Ml 



Z 3L c i 



M; 



T (v 2 -v v )-B- 



t'l 



I'll — + ^fijfkj^Vu k 



iB v Vi 



4 
.T 



■viv 2 - B, 



J kjJij 



viv Ui + - -^-vxv^ + ^fij(Xikj ~ X kij )v n v Vk + N'f, 



IB 



1 



■v 2 v Vi - B[, 



^(vl + vD-B^-B'^ 
v 2 



v 2 



4 

1 1 

(Mlh - \{g 2 - 9' 2 ){vl - vl + vl)5 l3 + 9 jV Ui v Vj + \r 3k n k v\ + n m 



(82) 
(83) 
(84) 
(85) 

(86) 
(87) 



+^fjiViV Vk (X k ii - X iH ) + -fuViv Uk {X kj i - Xjki), 



(89) 



1 

71 



+ ^Hk{X kj i - X jki )v 2 + ^ N ji V 2, 



(90) 



2^kifij v ^k v ^i + 2 



2"^ v y 1 2 

In Eqs. fl32|) and (0), the tadpoles conditions for the /i? and fields have been used: 



(91) 



■4 + /x 2 + i(s 2 + g' 2 ){vl - v\ + v 2 u )) Vl + { m + N'*)v Vi = 0, (92) 
(m 2 H2 + /i 2 + A - \{g 2 + g' 2 M - vl + vl))v 2 + Bv t + B[v Vi = 0. (93) 
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In Eqs. (|82| ) to ((?!]), v Ui are given by the corresponding tadpole conditions for the z/j fields: 

({MDji + fufi, + -(g 2 + g ,2 )(v 2 - v\ + vl)k^v u . + B[v 2 + N' i B v 1 + mVl = 0. (94) 

Without nonholomorphic terms, it is difficult to obtain a large cr — h + mixing without 
generating a large cl — h + mixing as well (which would induce an undesirably large W% — ► 
h + eL decay rate) or without requiring very fine tuning between the values of v v . and /ij. In 
the case of one family (putting all indices equal), this can be seen easily First, the A4 2 h -_ gc 
matrix element is proportional to the neutrino mass and hence very small. Second, the 
M 2 ,- ~ c matrix element, neglecting a small term proportional to the neutrino mass, is 



proportional to the Aij L _ gc matrix element. Hence it can be shown easily that it is not 
possible to have a sufficiently large M. 2 h ___ c matrix element (inducing £ of order 1CT 3 ) 
together with a sufficiently small — h~ mixing . The latter mixing gets a contribution 
~ £(yu/yU/)[yW| i _ g c/maa;(yW| i _g» , A^| c ,_g c )]. Now, in the case of three families, due to the 
A' e terms in M.\ L _^ C) both matrix elements are not any more proportional and the e l L — h + 
mixings can be made as small as necessary independently of the value of £. However, for 
values of §l and e c masses of the order 4 TeV or more (see section 4), a value of £ around 
10~ 3 requires that the A^^___ c matrix element is of order 10~ 3 x (4 TeV) 2 ~ (125 GeV) 2 
which implies very large values of \i\ and v n . Hence extreme fine tuning between the values 
of \ii and v Vi is needed to obtain a small enough neutrino mass in Eq. (Ill 
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Figure 1: Tree-level diagrams for (a) B' decay and (b) W% decay (through their B con- 
tent), and the one- loop (c) self-energy and (d) vertex diagrams for W' z decay which have 
absorptive parts of opposite lepton number. 
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Figure 2: Leptonic asymmetry as a function of z — M^, /T as obtained with the 
parameters given in the text, including all the contributions (solid); taking out the scatter- 
ing term (short-dashed); not considering in addition the inverse decay of the B' damping 
term (dotted); and without the inverse decay of the W% damping term (long-dashed). In 
the last case, since all damping terms have been taken out, the asymptotic result is just 
X L = en 7 /(2s). 
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Figure 3: Leptonic asymmetry Xi as a function of z = M^,JT as obtained with the 
parameters given by the set Eq. (pTFD, including all the contributions (solid); taking out 
the scattering term (short-dashed); not considering in addition the inverse decay of the 
B' damping term (dotted); and without the inverse decay of the W% damping term (long- 
dashed). In the last case, since all damping terms have been taken out, the asymptotic 
result is just Xl = en 1 /{2s). 
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Figure 4: Leptonic asymmetry Xi as a function of z = M^JT as obtained with the 
parameters given by the set Eq. (pRf), including all the contributions (solid); taking out 
the scattering term (short-dashed); not considering in addition the inverse decay of the 
B' damping term (dotted); and without the inverse decay of the damping term (long- 
dashed). In the last case, since all damping terms have been taken out, the asymptotic 
result is just Xl = en 7 /(2s). 
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Figure 5: One-loop diagram contributing to the sneutrino "Majorana" mass. 




Figure 6: One-loop diagram contributing to the neutrino mass, induced by the sneutrino 
"Majorana" mass. 
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